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Abstract. In this proceedings, we briefly review the methodology from [|T] for 
quantitatively constraining theoretical model parameters from experimental 
measurements including statistical and systematic uncertainties. We extend 
this methodology to additional parton energy-loss calculations for single in- 
clusive high pt particle suppression, and also extend the comparisons to di-jet 
observations. This is only the start of a process to give quantitative constraints 
on the quark-gluon plasma, and substantial theoretical uncertainties need to 
be reduced/resolved in a parallel path. 
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1. Introduction 

The field of relativistic heavy ion physics is undergoing a necessary transition from a 
field focused on the declaration of discovery - of thermalized partonic matter or the 
quark-gluon plasma, to one focused on the quantitative understanding of the unique 
properties of the medium created. This effort requires a more sophisticated and full 
treatment of experimental uncertainties, in particular systematic uncertainties and 
their underlying correlations. Additionally, the effort requires the reconciliation of 
variant theoretical approaches, or the discarding of some in favor of others. This 
proceedings focuses primarily on the first effort, though with some commentary on 
the second. 



2. Single Hadron High Suppression 

High pt tt°'s are suppressed in heavy ion reactions [[2], as shown in the left panel 
of Fig. [1] The experimental uncertainties are separated into three categories. Type 
A uncertainties are point-to-point uncorrelated and are shown as standard error 
lines. For this measurement, they are dominated by statistical uncertainties. Type 
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Fig. 1. (Color Online) The left panel shows the PHENIX ir° nuclear modification 
factor as a function of pr- The curves arc different calculations of the AMY+Hydro 
model corresponding to different input a s values. The right panel shows the corre- 
lation between the input a s coupling and the predicted nuclear modification factor. 

B uncertainties are point-to-point correlated and are shown as grey bars. They 
are dominated by the energy scale uncertainty in the Electro-magnetic calorimeters 
and have some contribution for px ~ 15 — 20 GeV/c from photon shower merging 
effects. Type C uncertainties are globally correlated (i.e. all points move by a 
common multiplicative factor). Note that this does not mean that all points move 
up and down together on a linear y-axis plot. This contribution is quoted as a ± 
12% global systematic uncertainty and has roughly equal contributions from the 
uncertainty in the nuclear thickness function Taa (as determined by a Glauber 
model) and from the proton-proton inelastic cross section absolute normalization. 
It is notable that these last two sources of uncertainty may prove quite difficult to 
reduce in the future. 

These Type B and C uncertainties correspond to one RMS deviations about 
the measured value. In the analysis that follows, these uncertainties are assumed 
to follow a Gaussian distribution with the corresponding RMS value. It is notable 
that if there are many contributions to the uncertainties, the Central Limit Theorem 
makes this assumption more solid. However, if the best fit results are a few standard 
deviations away in the systematic uncertainties, one should view the exact numerical 
fit result with some skepticism. Of course, if this is the case, the best fit typically is 
a very poor fit. It is also notable that if systematic uncertainties are quoted as 'full 
extent', then they cannot be added in quadrature and preserve their 'full extent' 
nature. 

Every publication of data in the field on which a full quantitative analysis is 
to be performed needs to explicitly quote these RMS uncertainty contributions in 
their appropriately labeled category. 
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Fig. 2. (Color Online) The modified \ 2 distribution as a function of the 
AMY+ Hydro model input parameter a s . 

The methodology is detailed in [[T], and here we just summarize the key result. 
One calculates a modified x 2 f° r a given data set with points {yi) and a given 
theoretical model with predictions [Xi as a function of parameter set p. One includes 
the possibility of systematic offsets given by the number of standard deviations for 
Type B (e b ) and for Type C (e c ): 



where &i is the uncertainty scaled by the multiplicative shift in yi such that the 
fractional error is unchanged under shifts 



Fundamentally one is determining if the penalty for moving the data by some num- 
ber of standard deviations in a systematic uncertainty is compensated by an overall 
reduction in the modified % 2 due to an improved statistical fit. 

An example comparison with a theoretical model is shown in Fig. Q] In this 
case, the /ij are the AMY+Hydro theoretical calculations [[3] as a function of the 
input parameters p, specifically being the coupling a s . In [ [3] , the authors utilize 
the AMY formalism for parton energy-loss and simulate the underlying medium 
with a hydrodynamic evolution model. The authors assert that "once temperature 
evolution is fixed by the initial conditions and evolution [by 3+1 dimensional hydro- 
dynamics], the coupling a s is the only quantity which is not uniquely determined." 
Shown in the left panel of Fig. [1] are calculations from this model for different input 
values of a s . The left panel shows how the nuclear modification factor at pt = 20 
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Model 


Constrained Medium Parameters 


p- value 


PQM 


(q) = 13.2 ^ 3 ,2 [1 stdev] and [2 stdev] GcV 2 /fm 


9.0% 


GLV 


{dN/dy) g i uon = U00tf 5 o [1 stdev] and [2 stdev] 


5.5% 


WHDG 


{dN/dy) gluon = 1400+^ [1 stdcv] and t<™ [2 stdev] 


1.3% 


zoww 


e = l-9to:s I 1 stdcv l and -ae I 2 stdev] GcV/fm 


7.8% 


AMY+Hydro 


a s = 0.280'lo:oia I 1 stdev] and [2 stdev] 


5.0% 


Linear 


b (intercept) = 0.168^ ^ [1 stdev] and [2 stdev] 

m (slope) = 0.0017±8;ggl| [1 stdev] and [2 stdev] (c/GeV) 


11.6% 



Table 1. Quantitative constraints on the medium parameters from various energy- 
loss models. 

GeV/c varies as a function of this input coupling value. We have then applied the 
full constraint method and show the modified y 2 as a function of a s in Fig. O 

The resulting constraint for the AMY + Hydro model and for a variety of parton 
energy-loss calculations and a linear functional fit are given in Table[T] It is critical 
to note that each constraint is assuming a perfect model calculation with only one 
unknown parameter, i.e. the uncertainty is from experimental sources only. 
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Fig. 3. (Color Online) For the PQM energy-loss model, the correlation between 
the (q)and the nuclear modification factor is shown. The left panel is a linear- linear 
display, and the right panel is a log-log display. 

There are a number of general observations that can be made at this point. 
In all of these calculations, there is a clear minimum in the modified x 2 and a 
reasonably well defined one- and two-standard deviation limit. In [0], the authors 
refer to the "fragility of high pt hadron spectra as a hard probe." The usage of this 
term "fragility" is ambiguous in the literature. Some simply refer to it with respect 
to the fact that as one increases the color-charge density or medium transport 
value (q), the nuclear modification factor Raa appears to saturate. Thus, one can 
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speculate that for some large value (for example (q)> 5 GeV 2 /fm), the uncertainty 
in determining Raa gets very large. However, viewing such trends on a linear y- 
axis scale can be quite deceptive. Shown in Fig. [3] are predictions from the Parton 
Quenching Model (PQM) [ [5] for different (q)values. In the right panel, shown 
on a log x-axis and log y-axis scale, is the modification factor versus (q). It is 
striking that it appears linear on this plot for (q)> 5. Thus, for a given fractional 
uncertainty in measuring Raa, one always gets the same fractional uncertainty on 
(q). Very similar results are obtained with other model calculations. 

If this is not the meaning of "fragility" (in a purely statistical constraint sense) , 
others refer to it in reference to the concept of surface emission bias. Imagine a beam 
of partons aimed at the corona of a dense medium. No measurement of the emitted 
particles provides information about the core of the medium (since the partons 
are not aimed there). However, if one had a model of the density distribution of 
the medium, then measuring (for example) the color-charge density of the corona, 
one can (via this distribution) learn about the density in the core. However, this 
knowledge is "fragile" in the sense that it depends more and more sensitively on 
the knowledge of the medium density distribution. It is interesting that in the 
paper discussing this "fragility" issue [ 0] , the authors employ a very unrealistic 
uniform cylindrical geometry, such that the density in the corona is identical to the 
density in the core. It is also notable that this geometry gets the distribution of 
hard scattering positions incorrect as well. 

3. Di-Jet Observables 




Fig. 4. (Color Online) The left panel shows the STAR dijet for per trigger away- 
side integrated yield (DdAulAv.Au) as a function of zt from d-Au and Au-Au. The 
right panel is the ratio (/,4a)- An additional 5% and 7% Type C uncertainty is 
including in the D^Au/AuAu and Iaa results, respectively [|8|. Also shown are the 
ZOWW theory calculations for different input energy-loss parameter (e ) values. 



6 



J.L. Nagle 




e GeV/fm 

Fig. 5. (Color Online) The modified \ 2 for the ZOWW calculations for the IauAu 
results (red line), the DauAu only (blue line), and the DauAu including the theo- 
retical scale uncertainty (magenta line). 

Regardless of the definition of "fragility" , it is always wise to search for addi- 
tional discriminating experimentally accessible observables. In [ 6J, they speculate 
that di-jet observables will be more sensitive to medium parameters. In the con- 
text of their calculation (ZOWW), they show that Iaa (the modification of per 
trigger away-side yields) has a steeper dependence on their energy-loss parameter 
Co (GeV/fm), than the single inclusive modification factor Raa- Thus, if one had 
identical experimental uncertainties on the two quantities, then Iaa would provide 
the stricter constraint. 

Shown in Fig. [4] are the experimental data from the STAR experiment [[7], and 
the calculations from the ZOWW model [[6|. If we apply a 7% Type C uncertainty 
to the STAR Iaa results [[8], the resulting constraint is eo = 2.9^ 'g [one std. dev.] 
and ig ? J [two std. dev.]. The '???' refers to the fact that within the parameter 
ranges available, there is no upper constraint. This is notably worse than the 
constraint from Raa given in Table [TJ In fact, the two results do not overlap each 
other within the one standard deviation uncertainties. The worse constraint from 
Iaa is due to the poor statistics in the d — Au reference (used instead of proton- 
proton) for the modification factor. This will be remedied in the near future with 
the much larger d — Au data set taken in 2007-2008. 

It is notable that in [ [6] , they show a rather tight eo constraint in a figure 
labeled Iaa- However, the actual constraint was derived from just comparing the 
Au-Au result alone (DauAu)- Thus, the constraint is only derived from the Au-Au 
data in the left panel of Fig. Q] The resulting modified x 2 's from Iaa and Da u Au 
alone are shown in Fig. [5l It is notable that the fit to the DauAu alone is a rather 
poor fit, which can be visually seen in Fig. |4] where the shape of the data points 
for zt > 0.4 do not match any of the theory curves. If one only fits DauAu, there 
is a definite additional theory uncertainty from the pQCD scale uncertainty in the 
NLO calculation (shown in [ 6 ). If this is included, the constraint is much looser, 
again shown in Fig. [5j Future high statistics data sets will allow for a more detailed 
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study. Also, the use of the variable zt may hide many details, and an optimal high 
statistics presentation should first include fine binning in trigger px for IauAu ■ 

4. Summary 

We summarize with a set of observations. On the experimental side, we have a 
well understood method for inclusion of statistical and systematic uncertainties. 
Experiments need to carefully quantify these Type A, B, and C uncertainties for 
all relevant measurements. It will be very interesting to see if larger p-p and d-Au 
data sets reconcile the Raa and Iaa constraints. 

On the theoretical side, one needs to resolve fundamental disconnects about 
whether perturbative calculations with modest a s fully describe parton energy loss. 
Currently the PQM results [ [5] indicate that the perturbative calculations fail. In 
either case, this must be reconciled at the appropriate scales with the picture of the 
bulk medium being strongly coupled (near-perfect fluid) and not being describable in 
a perturbative framework. Also, all calculations must include realistic geometries, 
fluctuations, and running of the coupling (if possible) so that the discussions and 
comparisons can focus on the more fundamental physics questions. 
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